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Abstract 

We apply a forward dispersion relation to the regeneration amplitude for kaon scattering 
on 12 C using all available data. The CPLEAR data at low energies allow the determination 
of the net contribution from the subthreshold region which turns out to be much smaller 
than earlier evaluations, solving a long standing puzzle. 
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1 Introduction 

Regeneration of neutral kaons is an excellent tool to study kaon-nuclear scattering [Jl]] 
and of crucial importance to experiments measuring CP violation in the kaon system [g]. In the 
latter context, CPLEAR [Q| has determined the kaon regeneration amplitude on C in the energy 
range between 0.56 and 0.9 GeV. In the present paper we explore the consequence of these data 
for strong-interaction physics analyzing the kaon-carbon scattering amplitude in the classical 
framework of forward dispersion relations (FDR) using in addition all the available information 
on KC scattering Q: the kaon regeneration on carbon total cross sections [JT3|— jTS]] , near 



forward differential cross sections JT7|], and the scattering length from kaonic atoms 
For high energies, we use the Regge model which is well supported experimentally [jOJ[T3|]. For 
the low energy application presented here this last aspect is not important. The energy range of 
the CPLEAR experiment is however well suited to obtain information on the strong-interaction 
physics occurring below the elastic threshold, often termed the unphysical region. Our results 
indicate that previous evaluations of FDRs have to be revised entirely in this region. Due to the 
new data situation, the net contribution from the unphysical region turns out to be much smaller 
and a simple coherent sum of elementary KN amplitudes is inadequate and even has the wrong 
sign as will be discussed below. 

2 Selection and treatment of data 

We define the regeneration amplitude as the difference between the nuclear amplitudes of 
forward scattering of K° and K° on C: 

/<-> H = /(K°C; u,9 = 0)- /(K°C; cu,6 = 0) (1) 

where cu is the total laboratory energy of the kaon. Note that this definition, which we choose 
for compatibility with earlier evaluations of FDR (an alternative definition with an extra factor 
1/2 in the r.h.s. of Eq. [j]is also used in the literature), has the opposite sign when compared to 
the A/ of our earlier regeneration analysis 0. 

Table 1 summarizes all experimental efforts to determine the regeneration amplitude in 
carbon. Most regeneration experiments determined modulus and phase of from de- 
cay rates behind a regenerator Two groups |H[|] a ^ so obtained a phase measure- 
ment from the time-dependent charge- asymmetry in semileptonic decays, thus independent of 
the CP violation phase 

Owing to the isospin symmetry of the 12 C nucleus, K 1 * 1 data may be used to extract infor- 
mation on the regeneration amplitude, if charge invariance is assumed: 

/H(o,) = /(K-C; u,9 = Q)- /(K+C; lu,6 = 0) (2) 

Then, with the use of the optical theorem, the difference between the total cross sections of K~ 
and K + is a measure of the imaginary part of f(~> [|I~5j-|r8|] (Table 2). One of these experiments 
[JTV|] also derived values for the real part from a measurement of the differential cross section 
dcr/dg 2 . 

Similarly, the K~C scattering length, a (K~C), extracted [ |20| ] from X-ray transition data 
obtained with kaonic carbon [ |19| ] may be regarded as a measurement of the K° scattering ampli- 
tude near zero kinetic energy. The corresponding amplitude for K°, having no inelastic channels 
open below threshold, must be purely real such that Im f^\m K a) = — Im ao(K~C). 

Where appropriate, the results have been adapted to the current world averages for the 



neutral kaon parameters, such as Am, rs and ?] + _, taken from Ref. pi]. Here it is important to 



1 ) All the data refer to natural carbon, with a 98.9% of 12 C content, to which the analysis applies. 
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note that some of these averages, in particular </>+_, are themselves influenced by regeneration 
measurements so that care must be taken when employing them to correct these very measure- 
ments. Such feedback effects, however, were found to be very small and may be neglected in 
our analysis. 

The corrected data are illustrated in Figs. [l|— 0, separately for measurements of modulus, 
phase, imaginary and real part of f(~>. Figure [3] reveals a conspicuous systematic disagreement 
between the total cross section results of Ref. [ |17| ] and the data of the other three experiments 



at the same energies. We exclude the data published in Ref. [17] from our fit, imaginary parts 



as well as real parts, as the latter directly depend on the former. For their impact see Section [4|. 

3 The dispersion relation and its parameterization 

In the framework of local relativistic field theories, analyticity of the amplitude has been 
proven to follow from local causality and it is therefore very natural to assume it as well for 
kaon scattering off nuclei. We use an unsubtracted dispersion relation as it is readily derived 
for antisymmetric amplitudes [|22|]. Since the precision of the data does not allow fitting several 
pole terms, they are represented by a single effective pole: 



2ur 2u; f°° , ,1m f^'[,,: 



Re /(">(„) = + — V I \> (3) 

where to denotes the total laboratory energy of the kaon and V stands for principal value. The 
effective pole term contains two parameters, position cu p and residue r. The lowest inelastic 
threshold lua = 183.6 MeV corresponds to A n X system at zero energy, see Fig. |5]a. The in- 
tegration over the imaginary part starts at cu^ and extends all the way to infinity (the physical 
region of K° and K° scattering on C corresponds to to > m K o). To carry out the integration, 
a parameterization of Im /(~) is needed over this range. Guided by the data and theoretical 
expectations, we have chosen the following parameterization: 

- The Regge model with one pole trajectory exchange is applied at energies above a 
certain energy uj k , giving 

Im/<->M =/3(^ ) Q sin(^[a + l]). (4) 

We choose the dimensional parameter cu = 1 MeV. With the slope parameter a ~ 0.42 
(see below), the rate of convergence of the dispersion integral @ at uj — ► oo is sufficient 
for reliable calculations in the near-threshold region. 

- In the intermediate energy range around 1 GeV, the total cross section data clearly indicate 
a resonance. We parameterize it for simplicity as a Gaussian on a background, the shape 
of which we describe by the same simple power law as employed at high energies, but 
with a different exponent: 



Im f(~\u) = b(u/u ) a sm(^[a + 1]) + Cexp j - 
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LU — U res 






a 





(5) 



The transition energy uj^ follows from a, (3, a and b if the function Im p)(u>) is to be 
continuous. In the fit, we will vary u; k in place of b as a free parameter. 
- The sensitivity of the CPLEAR data to the imaginary part of the regeneration amplitude 
is not sufficient for a parameter fit in the energy range below the first total cross section 
measurements in the 850 MeV region. Following Refs. p4] ] and p5|], we simply assume 



an essentially linear increase of Im p > from the physical threshold to 830 MeV total 
laboratory energy, in accordance with the data. 
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- At the physical threshold, we fix Im /(~) to the value corresponding to the measured 
K~C scattering length. The small error on this measurement is negligible in this analysis. 

- Virtually nothing is known about Im f(~>(u) = Im /(K°C; u, 9 = 0) below the elastic 
threshold, except that it must vanish at u^, the beginning of the cut. We find that any 
smooth and small contribution from the integral between u\ and rnxo in (0) is compati- 
ble with the experimental information, so that we may assume a linear behavior for this 
energy range as well. 

This gives a total of nine parameters to be determined in the fit. All of these parameters result 
from the interpolation of the data except for the effective pole, which represents the unphysical 
region below the elastic threshold. It turns out, however, that our evaluations are quite insensitive 
to the pole position cu p . We therefore fix it to the position of the A pole, given by the total energy 
needed to form a 12 C a resonance, 

s = rn 2 CA = m^o + m 2 -, + 2uj p m c , (6) 

which gives cu p = 172.3 MeV. 

Further parameters are needed for experimental reasons. As is evident from Fig. [1], there 
are systematic shifts in the normalization of the moduli measured by the high energy experi- 
ments (parameter /?) whereas there is good agreement on the slope parameter a. To prevent the 
systematic normalization uncertainties from affecting the result on a, we allow for three addi- 
tional correction factors, iV Car , AAib and N Roe for the Brookhaven [jlffi, Serpukhov [[TT|,^2|] and 



Fermilab [JT3p experiments, respectively, to be determined together with the other parameters 
in the fit. Since the later Fermilab publication [ [14| ] does not contain individual data points but 
only the value of a resulting from their analysis, we directly employ this value in our fit. 

4 Fit results and discussion 

We now vary the 1 1 parameters defined in the preceding section and compare at each 
data point the prediction for the real part from (BJ) with the experiment. Correlations between 
different measured quantities are taken into account where they have been reported, in particular 
the strong correlations between our own measurements of imaginary and real parts [Bp. Uncer- 
tainties in the kaon beam energy are neglected throughout. The weighted deviations are then 
summed up to a x 2 m the usual way. Minimization of this x 2 therefore reflects both the mea- 
sured imaginary and real parts of the amplitude and yields the following parameters (x 2 = 134.2 
for 127 degrees of freedom): 

r = -1.57 ±0.25 

C = (2.55 ± 0.18) fm 

u res = (1130 ±9) MeV 

a = (147 ±12) MeV 

uj k = (4570l|$) MeV (or b = (0.10 ± 0.02) fm) 

a = 0.61 ±0.03 

a = 0.424 ±0.005 

(3 = (0.35 ± 0.02) fm 

jV Car = 0.96 ±0.02 
iV A ib = 1.22 ±0.07 
A Roe = 1.07 ±0.03 
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The resulting fit functions are shown in Figs. [T]-Q, together with an optical model calculation 
[p6|]. The results in the near-threshold region are only weakly sensitive to the uncertainty in 
the high energy region. In particular, the total contribution of the Regge asymptotics (|]) from 
uj > cu k to the value Re f^\m K o) is about 0.9 fm, with the estimated error being less than 
0.1 fm. 

In order to discuss the physics below the elastic threshold we introduce the discrepancy 



function [22] 



AH = Re - Re (a,) (7) 

where Re f^ai 00 ) represents the dispersion integral (Eq. from the physical regionF] (u > 
m K o). If the unphysical region is well represented in the physical region by the effective pole 
assumed in Eq. |]then the quantity 

D(u) = — ^ (8) 

should be independent of the energy and equal to the effective residue^], 

D(u) = r. (9) 

In Fig. 0we show D(u) for to p from Eq. || corresponding to a 12 Ca system. The function D(u) is 
not sensitive to the exact pole position and Fig. shows no energy dependence which confirms 
the validity of the effective pole ansatz. The effective residue turns out to be small and negative 
while previous attempts to determine it from FDR led to large and positive values: 

-1.57 ±0.25 a) this analysis 

23.2 ±5.8 b) Ref. [24] (10) 
12.2 ± 3.0 c) Ref. 



The main reason for this is that both of these earlier analyses completely relied on the data of 
Ref. JT7| ] which we have discarded in view of their incompatibility with all other experiments. 
Given the large error bars of these data, their inclusion into our fit hardly changes the pole 
residue but increases the x 2 to 197.6 for now 133 degrees of freedom. Without the CPLEAR 
data for the real part, we would get a large positive residue of order of r = 12, as can be seen 
from Fig. consistent with the earlier analysis [|25l|. 



Our main result therefore is that the residue is small. In our analysis even its negativity 
is stable against considerable variations of a smooth background in the unphysical region and 
against variations of the imaginary part in the physical region below 830 MeV. However, since 
the contribution of the effective pole is not large where the real parts are measured by CPLEAR, 
0.56 < uj < 0.9 GeV, the sign is less certain than the magnitude. 

From first principles the sign of the residue is directly related to the parity of the system 
of the particles exchanged $23j\. For the elementary K ± p scattering, hyperon exchange leads to 
a positive residue (similar to 7rN scattering). Hence if the unphysical region for K 12 C scattering 
were dominated by the superposition of elementary hyperon exchange of Fig. 0a, a positive 
residue would be expected. If, however, an additional pion in a relative S-wave is exchanged 
(Fig. 0b) or if the hyperon in Fig. 0a is in P-wave relative to the spectator nucleons, then the 

2 ) To avoid a singularity at u) = m^o in fp h y S (^) we use a linear parameterization of Im (uj) for u < m^o 
as described in Section [3| Its contribution to A(u>) is negligible in the energy range of CPLEAR data. 



3 ) More elaborated methods (see [ 28 1) do not help in our case because of the scarcity of data. 
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opposite parity results and the residue is negative. Our analysis clearly indicates that there are 
large cancellations between these contributions with different parities since the residue is stably 
small. The smallness of the effective residue is actually very gratifying since the earlier results 
[ p4| , |25| ] claiming a large and positive residue are very hard to understand. As was already noticed 
by Dumbrajs [ |25ll , the impulse approximation has no chance to work. He observed that the 
measured imaginary part (total cross section) at low energies are much smaller than the impulse 
approximation while the result for the effective residue (|T0|c) is about 50% larger than the sum 
of the elementary residues for A and E exchanges. The value (|10"|b) from p3| ] is even bigger. 
Our result (|T0|a) therefore solves a long standing puzzle. 

The threshold region deserves some comment. In view of the large and repulsive K~C 



scattering length [g0| we have also used the near threshold parameterization of [[240 where 
this feature is related to a narrow peak in the imaginary part below threshold, see Fig. 3. 
This local structure near threshold does not affect our parameters for the effective pole, as 
we have checked. Our evaluation gives for the real part of the regeneration amplitude at thresh- 
old Re/(~)(m K o) = 2.9 fm. Using the real part of the K^C scattering length Rea (K~C) = 
2.9 fm from Ref. p0| ] we obtain the real part of the K + C scattering length Reao(K + C) = 
Re f^~\m K o) + Rea (K"C) = 5.8 fm. The value of a (K + C) should be considered as an esti- 
mate: it is sensitive to the parameterization of Im f(~> (u) near the elastic threshold because of a 
steep energy dependence of Re (to) (see Fig.|j). For comparison, using an optical model, we 
obtain a (K + C) ~ 1.5 fm. Obviously there is a partial cancellation between the contributions 
from K~ and K + for /(~) since both scattering lengths are repulsive. While these findings are 
interesting for the discussion of scattering lengths and the detailed behavior of the amplitudes 
near threshold, these local structures are irrelevant for our discussion of the effective pole which 
is dominated by much larger energy scales (the energy range of the CPLEAR experiment). 

In summary, we conclude that the combined information on the K 12 C and K 12 C ampli- 
tudes is very well described by a forward dispersion relation containing only a small effective 
pole contribution from the subthreshold region, consistent with theoretical expectations. Earlier 
evaluations claiming a large contribution from the unphysical region can be excluded on the ba- 
sis of the measured regeneration amplitude from CPLEAR in combination with the information 
on this reaction from all other consistent experimental data. These results are of relevance in the 
analysis of nuclear scattering and of CP- violation experiments at low energy [|29|]. 
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Table 1: Carbon regeneration experiments. 
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Likelihood plots in the complex plane are given. 

b Published data are analyzed using the value of 0_| resulting from the compilation of Ref. [^T|. 

6 Only a fitted value for the power law exponent a is published. 



Table 2: Carbon Pv^ experiments. 

group Ref. p K measured derived x 2 /dof 

[GeV/c] quantities part of p' in fit 



Bugg et al. 
Abrams et al. 
Gobbi et al. 

Afonasyev et al. 
Backenstoss et al. / 
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■ Carithers et al. [10] 
□ Bott-Bodenhausen et al. [7] 
o Bohm et al. [5] 
Christenson et al. [4] 
Angelopoulos et al. [3] 
combined with Bugg et al. [15] 
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Figure 1 : Modulus of the regeneration amplitude /( experimental results in comparison with 



our fit (solid line). The data points of Schwingenheuer et al. are taken from Ref. 027| ] and are 
shown just for visualization (only their published value for a enters the fit). Note the double- 
logarithmic scale. 
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Roehrig et al. [13] 
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Figure 2: Phase of the regeneration amplitude p '•. experimental results in comparison with 



our fit (solid line). The data points of Schwingenheuer et al. are taken from Ref. p7Q and are 
shown just for visualization (only their published value for a enters the fit). The data of Refs. [j^] 
and [0] as well as those of Refs. [|]] and [§] have been averaged for clarity. Note the logarithmic 
scale. 
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co(K u ) [GeV] 



Figure 3: Imaginary part of the regeneration amplitude f(~>: experimental results in comparison 
with our fit (solid line) and the optical model calculation of Ref. P^| ] (dashed line). The dotted 
line represents the sub-threshold parameterization employed in Ref. []24|]. The data of Ref. [ |T7| ] 
are discarded in the fit. Note the linear scale. 



10 




co(K u ) [GeV] 

Figure 4: Real part of the regeneration amplitude /' experimental results in comparison with 
our fit including the pole at u p = 172.3 MeV (solid line) and the optical model calculation of 
Ref. p6| ] (dashed line). The data of Ref. [ |T7| ] are discarded in the fit. Note the linear scale. 
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Figure 6: The function D(u) as derived from our data (full circles) and from the data of 
Ref. n (see text). 
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